EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP/2012-244 
2012/09/05 



CMS-EXO-1 1-099 



Search for pair produced fourth-generation up-type quarks 
in pp collisions at = 7 TeV with a lepton in the final state 



The CMS CoUaboratioif] 



Abstract 



The results of a search for the pair production of a fourth-generation up-type quark 
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momentum, and at least four hadronic jets, of which at least one must be identified as 
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is observed. Upper limits for the t'F production cross section at 95% confidence level 
are set as a function of t' mass, and t'-quark production for masses below 570 GeV is 
excluded. 
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1 Introduction 

The standard model (SM) of particle physics includes three generations of fermions How- 
ever, the existence of a fourth generation is not excluded by precision electro weak data IH. 
Mass scales above 550 GeV are interesting because the coupling of a fourth-generation up-type 
quark (referred to as t' in this paper) to the Higgs field becomes large, weak interactions start 
to become comparable to strong interactions, and perturbative calculations begin to fail [5J . In 
addition, some models of physics beyond the SM HI IZl predict the existence of top-like quarks 
that cancel the divergent corrections of t-quark loops to the Higgs boson mass. Such quarks 
give rise to the same final-state signature described below for fourth-generation t' quarks. 

The mass splitting between the t' quark and the corresponding fourth generation down-type b' 
quark is favoured to be smaller than the mass of the W boson ISl- ITOt . In this case, the t' quark 
cannot decay to Wb'. Assuming that the pattern of quark mixing observed in the CKM matrix 
extends to the fourth generation, the dominant t'-quark decay mode would be t' — > Wb, and 
the lifetime of the t' would be short in an experimental sense. 

We present the results of a search for the strong production of a t'F quark pairs, with t' decaying 
into W+b and F to W^b, in proton-proton collisions at = 7TeV using the Compact Muon 
Solenoid (CMS) detector. The search strategy requires that one of the W bosons decays to 
leptons (ev or jiv) and the other to a quark-antiquark pair. The branching fraction into this 
final state is about 15% for each lepton flavour. We select events with a single charged lepton, 
missing transverse momentum, and at least four jets with high transverse momenta (px)- 

Previous searches for t' quarks in this final state give lower limits for the mass of the t' quark 
of 358 GeV llllp] at the Tevatron and 404 GeV ifil at the Large Hadron Collider (LHC). 

There are SM processes that give rise to the same signature, most notably tt and W+jets pro- 
duction. The present search considers a t' quark with a mass larger than the SM t quark. We 
utilize two variables to distinguish between signal and background. The first is Hj, defined 
as the scalar sum of the transverse momenta of the lepton, the missing transverse momentum, 
and the four jets from the decay of the t' and F quarks. The second variable is the t'-quark mass 
Mfit, obtained from a kinematic fit of each event to the process t'F — ?► WbWb — ?► £i/bqq'b . We 
use the two-dimensional distribution of Hj versus Mgt to test for the presence of a signal for 
t'F production in the data. 

We categorize events according to the flavour of the lepton. Events with an identified electron 
(muon) are classified as e+jets (/i+jets) events. The analysis procedures for the two channels are 
kept as similar as possible, with small differences mainly driven by the different trigger condi- 
tions. Finally, a combined statistical analysis of both channels is performed and upper limits 
for the t't' pair production cross section and a lower limit on the t'-quark mass are derived. 

2 CMS detector and data samples 

The CMS experiment uses the following coordinate system. The z axis coincides with the axis 
of symmetry of the detector, and is oriented in the anticlockwise proton beam direction. The 
X axis points towards the center of the LHC ring and the y axis points up. The polar angle 
9 is defined with respect to the positive z axis, and (p is the azimuthal angle. The transverse 
momentum of a particle or jet is defined as its momentum times sin 6, and pseudorapidity is 
}] = -ln[tan(0/2)]. 

The characteristic feature of the CMS detector is the superconducting solenoid, 6 m in diameter 



2 



3 Event reconstruction 



and 13 m in length, which provides an axial magnetic field of 3.8 T. Inside the solenoid are a 
multi-layered silicon pixel and strip tracker covering \ri\ < 2.5 to measure the trajectories of 
charged particles, an electromagnetic calorimeter (ECAL) made of lead tungstate crystals and 
covering \ f]\ < 3.0, a preshower detector covering 1.65 < |f/| < 2.6 to measure electrons and 
photons, and a hadronic calorimeter (HCAL) made of brass and scintillators covering \f]\ < 3.0 
to measure jets. Muons are identified using gas-ionization detectors embedded in the steel 
return yoke of the solenoid and covering < 2.4. Extensive forward calorimetry comple- 
ments the coverage provided by the barrel and endcap detectors. The CMS detector is nearly 
hermetic, allowing the measurement of the transverse momentum carried by undetected par- 
ticles. A detailed description of the CMS detector can be found elsewhere I|l4) . 

We use data collected in 2011, corresponding to an integrated luminosity of 5.0 fb~^ for the 
e+jets channel and 4.9 fb^^ for the /^+jets channel. The triggers for the e+jets data required at 
least one electron candidate with a px threshold that varied between 25 and 32 GeV according 
to the average instantaneous luminosity. When the LHC instantaneous luminosity increased, 
three central jets with pj > 30 GeV and \ f]\ < 2.6 were also required. The triggers for the f/+jets 
channel required at least one muon candidate with a pj threshold that varied between 30 and 
40 GeV. No requirements were made on jets in the triggers for the /^+jets events. 

We model the t'F signal and SM background processes using Monte Carlo (MC) simulations. 
The t'F signal events are generated for t' masses from 400 to 625 GeV in 25 GeV steps. The 
following SM background processes are simulated: tf production; single t-quark production 
via the tW, s-channel, and f-channel processes; single- and double-boson production (W+jets, 
Z+jets, WW, WZ, and ZZ). All of these processes, except the dominant tt production, are 
collectively referred to as electroweak (EW) background. 

The single t-quark production is simulated with the POWHEG event generator |[T5l - [T7l . All 
other processes are simulated with the madevent/MadGraph [18] programs. The PYTHIA 
program |lT9t is then used to simulate additional radiation and the fragmentation and hadroni- 
zation of the quarks and gluons into jets. The generated events are processed through the CMS 
detector simulation based on Geant4 [20J. Up to 20 minimum-bias events, generated with 
PYTHIA, are superimposed on the hard-scattering events to simulate multiple pp interactions 
within the same beam crossing. The MC events are weighted to reproduce the distribution of 
the number of vertices per event in the data (the average number of vertices per event is 8). 

The simulated samples for the t'F signal correspond to an integrated luminosity of between 100 
and 2500 fb^^ for each value of t'-quark mass. Samples for the background processes giving the 
largest contributions correspond to 22 fb^^ for the tt sample and 2.5 fb^^ for W+jets. 

3 Event reconstruction 

Events are reconstructed using a particle-flow algorithm ['2T'-'231. The particle-flow event re- 
construction consists in reconstructing and identifying each single particle with an optimized 
combination of all subdetector information: charged tracks in the tracker and energy deposits 
in the ECAL and HCAL, as well as signals in the the muon system and the preshower detector. 
This procedure categorizes all particles into five types: muons, electrons, photons, charged and 
neutral hadrons. The energy calibration is performed separately for each particle type. 

Electron candidates are reconstructed from clusters of energy deposited in the ECAL. The clus- 
ters are first matched to track seeds in the pixel detector. The trajectory of the electron can- 
didate is reconstructed using a dedicated modelling of the electron energy loss. Finally, the 
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particle-flow algorithm further distinguishes electrons from charged pions using a multivari- 
ate approach p3| . 

Muon candidates are identified by reconstruction algorithms using signals in the silicon tracker 
and muon system. The tracker muon algorithm starts from tracks found in the tracker and then 
associates them with matching signals in the muon chambers. The global muon algorithm 
starts from standalone muons and then performs a global fit combining signals in the tracker 
and muon system. 

Jets from the fragmentation of quarks and gluons are reconstructed from all particles found 
by the particle-flow algorithm using the anti-A;x jet clustering method [24] with the distance 
parameter of R = 0.5, as implemented in Fastjet version 2.4 [25r.27,|. Small corrections |28l 
are applied as a function of tj and px to the reconstructed jet energies. 

A jet is identified as originating from a b quark using the combined secondary vertex (CSV) 
algorithm Il29ll , which provides optimal b-tagging performance. The algorithm is based upon a 
likelihood test that combines information about the impact parameter significance, secondary- 
vertex reconstruction, and jet kinematics. The small differences in the performance of the b- 
tagging algorithm in data and MC simulation are accounted for by data/MC scale factors. 
This is done by randomly removing or adding b tags on a jet-by-jet basis, using the and 
(/-dependent scale factors discussed in P29| . 

The missing transverse momentum in an event is defined as the negative vector sum of the 
transverse momenta of all objects found from the particle-flow algorithm. 

The vertex with the highest sum of p\ of all associated tracks is taken as the primary vertex of 
the hard collision. 

4 Event selection, signal and background estimation 

For this analysis we use an event selection similar to that adopted previously for tf events in 
the lepton+jets channel li30J . To reduce the background from tf production, we apply higher jet 
px thresholds. 

Charged leptons from W-boson decays, which are themselves originating from decays of heavy 
t-quark-like objects, are expected to be isolated from nearby jets. A lepton isolation variable is 
calculated by summing the transverse momenta of all reconstructed particle inside a cone de- 
fined as AR = ^J {^<pY + {^vY — 0-3, where and A// are the azimuthal angle and pseudo- 
rapidity differences with respect to the lepton direction. The lepton isolation variable is equal 
to this sum divided by the lepton's px- 

Events with exactly one isolated lepton and at least four jets with < 2.4 are selected. Jets 
that are within a cone of AR = 0.3 around the lepton direction are not considered. At least 
one jet must be identified as originating from a b quark. The thresholds for the lepton px 
are driven by the trigger requirements described in Section |2j The lepton track must have an 
impact parameter transverse to the beam direction with respect to the primary vertex of less 
than 0.02 cm and along the beam direction of less than 1 cm. The missing px in the event must 
be greater than 20 GeV. 

The selection of the e+jets events requires exactly one electron with px > 35 GeV and \y]\ < 2.5, 
electron isolation < 0.1, and at least four jets with px > 120, 90, 50, and 35 GeV. The selection 
for the /^+jets channel requires exactly one muon with px > 35 GeV or px > 42 GeV for two 
running periods with different trigger conditions, \rj\ < 2.1, muon isolation < 0.125, and at 
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least four jets with pj > 120, 90, 30, and 30 GeV. The thresholds for the two highest-px jets are 
selected to maximize the signal-to-background ratio. The thresholds for the lepton pj and the 
third and fourth highest-pj jets are driven by the trigger conditions. 

Table [l] lists the number of observed and expected events for the various background sources 
after selection. The expected numbers of background events are calculated from the cross sec- 
tions and integrated luminosities given in the table. The cross section for tt production is taken 
from a previous CMS measurement II30I . All other cross sections are computed with the MCFM 
program [31 1. In the case of the e+jets channel, the small multijet background is estimated 
from data by fitting the missing-pj distribution with shapes predicted by the MC simulation. 
The uncertainties shown include systematic uncertainties in the efficiency and acceptance. Un- 
certainties are strongly correlated for all sources. Uncertainties in the cross sections and the 
integrated luminosity are not included. 

The fraction of tt events retained by our selection is 0.7% for the /i+jets channel and 0.5% for 
the e+jets channel. 

Table 1: Background cross sections, number of events observed and background events pre- 
dicted for the e+jets and /^+jets samples. The predicted numbers of events are normalized to 
the integrated luminosity (except for the multijet events in the e+jets channel, see text). 







e+jets 




Integrated luminosity 




4.98 fb"^ 


4.90 fb"^ 


Background process 


Cross section 


Events 


Events 


tt 


154 pb 


3950 ± 490 


5460 ± 670 


W+jets 


31 nb 


462 ± 55 


750 ±110 


Single-t production 


85 pb 


208 ± 24 


336 ± 45 


Z+jets, WW, WZ, ZZ 


3.1 nb 


49 ±8 


69 ±11 


Multijets 




78 ±9 


5±5 


Total background 




4750 ± 560 


6620 ± 800 


Total observed 




4734 


6448 



Table |2] shows the theoretical cross sections for the signal process for various t'-quark masses, 
along with the efficiencies of the event selection for the e+jets and f/+jets channels and the 
expected numbers of signal events. The t'F production cross sections are computed using 
HATHOR [32 J. The efficiencies include the branching fraction of the t'F system into a single- 
lepton final state, which can be obtained from the branching fractions for W — )■ £v and W — ?► 
qq'. The uncertainties quoted are the statistical uncertainties from the MC simulation. 

5 Mass reconstruction 

We perform a kinematic fit of each event to the t'F — )■ WbWb — ?► £i/bqq'b process. There are two 
steps in the reconstruction of the t'-quark mass: the assignment of reconstructed objects to the 
quarks, and a kinematic fit to improve the resolution of the reconstructed mass of the t'-quark 
candidates. The four-momenta resulting from the kinematic fit of the particles in the final state 
must satisfy the following three constraints, where m is the mass of the corresponding particles, 
Mw is the W-boson mass. Mat is a free parameter in the fit (reconstructed t' mass), and £ stands 
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Table 2: Theoretical cross sections II32I , selection efficiencies, and numbers of expected events 
for the t'F signal with different t' masses in the e+jets and /^+jets channels. The efficiencies 
include the branching fraction of the t'F system into a single-lepton final state. 



Mf (GeV) 


Cross section (pb) 


e+jets eff. (%) 


Events 


t/+jets eff. (%) 


Events 


400 


1.41 


4.3 ±0.1 


302 


5.4 ±0.1 


373 


425 


0.96 


4.4 ±0.1 


210 


5.6 ±0.1 


263 


450 


0.66 


4.7 ±0.1 


155 


6.0 ±0.1 


194 


475 


0.46 


4.7 ±0.1 


108 


6.1 ±0.1 


137 


500 


0.33 


4.8 ±0.1 


79 


6.2 ±0.1 


100 


525 


0.24 


4.7 ±0.1 


56 


6.4 ±0.1 


75 


550 


0.17 


4.9 ±0.1 


41 


6.5 ±0.1 


54 


575 


0.13 


4.7 ±0.1 


30 


6.6 ±0.1 


42 


600 


0.092 


4.7 ±0.1 


22 


6.6 ±0.1 


30 


625 


0.069 


4.8 ±0.1 


16 


6.5 ±0.1 


22 



for electron or muon: 

m{£v) = Mw, (1) 

m(qq') = Mw, (2) 

m{ivh) = m(qq'b) = Mfit. (3) 

Here £, v, b denote either particle or antiparticle. 

The reconstructed objects in the event are the charged lepton, the neutrino, and four or more 
jets. For the neutrino, only its transverse momentum can be measured as the missing trans- 
verse momentum in the event. The z component of the neutrino momentum can be deter- 
mined with two solutions from the kinematic constraints. The four quarks in the final state 
manifest themselves as jets and their momenta are measured. Thus, all but one of the momen- 
tum components of the considered final system are measured. With one unknown and three 
constraints, a kinematic fit is performed by minimizing the ^(^ computed from the difference 
between the measured momentum of each reconstructed object and its fitted value, divided by 
its uncertainty. 

We have studied different strategies for pairing the observed jets with the four quarks from 
the decay of the t' and F quarks to find the best separation between the t'F signal and the tt 
background. In events with exactly four jets, we consider all possible jet-quark assignments. 
To reduce the number of combinations, we choose only those in which at least one b-tagged jet 
is assigned to a b quark from the t'F decay. In events with more than four jets, we take the five 
jets having the highest px values, and consider all combinations of four out of these five jets. 
In each event, the kinematic fit is carried out for each jet-quark assignment, and the jet-quark 
assignment with the smallest value is chosen. This procedure selects the correct jet-quark 
assignment in 36-40% of the simulated t'F events over a t'-quark mass range of 400-625 GeV for 
all jet multiplicities together. For tt events this fraction is much lower, about 19%, because the 
jets from the decays of the t and t quarks are softer than from t' and F decays and, therefore, are 
less likely to be among the five highest-px jets in the event. The value does not distinguish 
the t'F signal from the tt background because both processes satisfy the fit hypothesis, but 
using the smallest value for each event does increase the fraction of correct quark assignments. 
Since a restriction on does not improve the signal-to-background ratio, no such restriction is 
applied. 
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Figures [T] and |2] show the two-dimensional Hj versus Mfit distributions for the data, tt sim- 
ulation, the other simulated backgrounds, and the t'F simulation with a particular t' mass of 
550 GeV in the e+jets and //+jets channels, respectively. Figure |3] shows the corresponding Mm 
and Hj projections. The integrated luminosities given in Table [T] are used for the normaliza- 
tion of the background processes. The data are found to be in agreement with the predicted 
background Mm and Hj distributions. The tf events that pass the selection criteria either have 
high-px t and t quarks that produce high-px jets in their decays or they have high-px jets from 
initial-state gluon radiation. The former class of events is responsible for the relatively nar- 
row peak in the Mm distribution at the t-quark mass. The Mm distribution of the latter class 
of events is broad and typically populates the region above the t-quark mass, leading to the 
observed high-mass tail in the Mfit distribution. 

6 Computation of tT cross section limits 

The two-dimensional distributions of Hj versus Mm, such as those shown in Figs. [T] and |2| are 
used to search for a t'F signal in the data. Simulated t't' signal distributions are produced for 
t' masses from 400 to 625 GeV in 25 GeV steps. We do not use the two-dimensional histograms 
directly because it is not possible to simulate enough events to adequately populate all bins of 
the distributions for both signal and background. Therefore, we employ a new procedure that 
combines bins. 

All the background distributions are added together to obtain the expected background event 
yield in each bin of the Hj versus Mfit histogram. Then the projections of the two-dimensional 
signal and background histograms onto the Hj and Mfit axes are separately fitted with ana- 
lytic functions. Next, we compute the expected signal-to-background (s/b) ratio for each two- 
dimensional bin as the product of the values of the two one-dimensional-bin fit functions for 
the signal and for the background at the bin center. All two-dimensional bins are then sorted 
in increasing order of the expected s/b ratio, which we call s/b rank. Using the fit functions 
for calculating the s/b ordering rank removes biases due to statistical fluctuations in the simu- 
lated signal and background samples. These functions are used only to define the ordering of 
the bins. 

We then merge the two-dimensional bins that are adjacent after ordering by s/b ratio so that 
the fractional statistical uncertainty of both the signal and the background predictions is below 
20% in all bins. We select the 20% value as a compromise between two tendencies. Values above 
20% would give better signal sensitivity, but they could lead to significant biases in the t'F cross 
section determination. Figure |4] shows the colour-coded maps of the merged bins obtained for 
the simulation of a t' quark with a mass of 550 GeV. The colour represents the rank of the bin 
in the s/b ordering. A higher rank corresponds to a higher s/b value. 

In Fig. |5| the number of events in the merged bins is plotted versus s/b rank. In these his- 
tograms, signal events will predominantly cluster towards the right, and background events 
towards the left. These one-dimensional histograms are used as input to the t'F cross section 
computation, and we will refer to these distributions as templates in the following. The data 
agree with the predicted background distributions in Fig. [sj with no evidence for a t' signal. 
Thus, we use the results to set upper limits on the t'F cross section as a function of t' mass. 

The computation of the limits for the t'F cross section uses the CLg criterion [3311341. The first 
step is to perform a likelihood fit to the data. We group the background in two components: 
the larger one due to tt production and the smaller one from all EW processes (W+jets, Z+jets, 
single-t, and diboson production) and from multijet processes. Each background component is 
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Figure 1: Hj versus Mgt for the e+jets channel from data (top left), and simulations of tt pro- 
duction (top right), other backgrounds (bottom left), and t'F production (bottom right) for 
Mf = 550 GeV. 
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Figure 2: Hj versus Mgt for the /i+jets channel from data (top left), and simulations of tt pro- 
duction (top right), other backgrounds (bottom left), and t'F production (bottom right) for 
Mf = 550 GeV. 



8 



6 Computation of t't' cross section limits 



CMS is = 7TeV L = 5.0fb^ CMS /s = 7TeV L = 5.0fb^ 




200 400 600 800 500 1000 1500 

Mfjt [GeV] [GeV] 



CMS is = 7TeV L = 4.9fb^ CMS /s = 7TeV L = 4.9fb^ 




200 400 600 800 500 1000 1500 

Mfit [GeV] [GeV] 

Figure 3: Distributions of Mfit (left) and Hj (right) for the e+jets (top) and ^+jets (bottom) 
channels. The data are shown as points, the simulated backgrounds as shaded histograms, and 
the expected signal for a t' mass of 550 GeV as dashed histograms (multiplied by a factor of 50 
to improve visibility). 
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Figure 4: Map of the merged bins in the Hj versus Mm plane for a t' quark with a mass of 
550 GeV for the e+jets (left) and f/+jets (right) channels. The colour represents bins merged 
according to increasing signal-to-background (s/b) ratio. The vertical colour axis is labelled by 
s/b rank and corresponds to the bin index of the one-dimensional histograms of Fig.|5j 
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Figure 5: Number of events per bin in the two-dimensional Hj versus Mgt histogram after bin 
merging, as a function of the signal-to-background (s/b) rank for the e+jets (left) and f^+jets 
(right) channels. The data are shown by the points, the simulated tt and other background 
distributions by the histograms, and the prediction for a t'F signal with a t' mass of 550 GeV by 
the dotted lines (multiplied by a factor of 50 for easier viewing). 



normalized to its expected yield and multiplied with a scale factor that is a free parameter in 
the fit. The t'F cross section, a, is also a free parameter in the fit. The following likelihood ratio 
is used as the test statistic: 

1 1 if a < a. 

Here, L{q\cr,<x) is the likelihood of the data having the value q for the parameter of interest 
and the nuisance parameters a. The nuisance parameters account for effects that give rise to 
systematic uncertainties in the templates and include the normalizations of the background 
components. The likelihood reaches its maximum when a = a and a = a. The symbol occr 
refers to the values of the nuisance parameters cc that maximize the conditional likelihood at a 
given value of cr. 

Using the asymptotic approximation for the test statistic described in ||35| , the probability to ob- 
serve a value t for the likelihood ratio that is larger than the observed value tgjjs is determined. 
This is done by producing samples of "pseudo-experiments" in which the expected numbers of 
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signal and background events are allowed to vary according to their statistical and systematic 
uncertainties. For the pseudo-experiments generated with background only, this probability is 
denoted by CL},. For pseudo-experiments with a cross section a for the t'F signal, this probabil- 
ity is denoted by CLg+j, {a), which is a function of a. The upper limit at the 95% confidence level 
(CL) for the t'F cross section is the value of a for which CLg = CLg^},/ CL^ = 0.05. To determine 
the limits for both lepton channels combined, we simultaneously fit the histograms from both 
channels, accounting for correlations among the nuisance parameters, and then apply the CL, 
method described above. 

7 Systematic uncertainties 

The signal and background predictions are subject to systematic uncertainties. Below, we de- 
scribe all sources of systematic uncertainties that have been considered. They can be divided 
into two categories: uncertainties that only impact the normalization of the signal and back- 
ground templates, and uncertainties that also affect the shape of the distributions. 

The uncertainties in the tt cross section, electroweak and multijet background normalizations, 
integrated luminosity, lepton efficiencies, and data/MC scale factors affect only the normaliza- 
tion. 

The uncertainty on the cross section for tt production is taken from the CMS measurement 
of 154 ± 18pb at = 7TeV ||30|. The predicted yields of the EW and multijet backgrounds 
are determined as described in Section |4] A 50% uncertainty is assigned to the sum of these 
two backgrounds in the likelihood fit to the data in order to account for the uncertainty in the 
acceptance and the W+jets normalization. 

The integrated luminosity affects the normalization of the t'F signal and the background tem- 
plates in a correlated way. The integrated luminosity is known to a precision of 2.2% [36 1. 

Trigger efficiencies, lepton identification efficiencies, and data/MC scale factors are obtained 
from data using decays of Z bosons to dileptons. Their uncertainties are included in the se- 
lection efficiency uncertainty. They amount to 2% for the /^+jets channel and 3% for the e+jets 
channel. 

Uncertainties that affect the shape of the distributions include those on the jet energy scale, jet 
energy resolution, missing-px resolution, b-tagging efficiency, number of multiple pp interac- 
tions, factorization/ renormalization scale Q, matrix-element/ parton-shower matching thresh- 
old Il37ll , and initial- and final-state radiation. To model these uncertainties, we produce addi- 
tional templates by varying the nuisance parameter that characterizes the systematic effect by 
±1 standard deviation. To determine the signal and background templates used in the fit for 
any value of the nuisance parameter, we interpolate the content of each bin between the varied 
and nominal templates. This procedure is often referred to as vertical morphing. 

The energy of all jets is obtained using the standard CMS jet energy calibration constants Il28| . 
The sum of the four-momenta of the jets is 100% correlated with the measured missing px- 
The jet energy scale uncertainty affects the normalization and the shape of the Hj vs. Mgt 
distribution. This is taken into account by generating Hj vs. M^t distributions for values of the 
jet energy scaled by ±1 standard deviation of the rj- and px"dependent uncertainties from Il28l . 

The energy resolution of jets in the simulation is better than in the data. Therefore, random 
noise is added to the jet energies in the simulation to worsen the resolution by 10%, to match 
the actual resolution of the detector. To estimate the corresponding uncertainty, the analysis 
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is performed without smearing and with 20% smearing. The missing-px resolution is also 
simultaneously corrected for this effect. 

The systematic uncertainty from the b-tagging efficiency is estimated by varying this efficiency 
by ±1 standard deviation taken from [29 J. 

To evaluate the uncertainties related to the modelling of multiple interactions in the same beam 
crossing, the average number of interactions in the simulation is varied by ±8% relative to the 
nominal value. 

The uncertainty in the factorization/ normalization scale Q, used for the strong coupling con- 
stant tx.s{Q^), is estimated by using two sets of simulated tt samples in which the Q value is 
increased and decreased by factors of two relative to the nominal value. 

The uncertainty arising from the threshold for matching between matrix elements and parton 
showers ['37] is estimated using two simulated tt samples generated with the matching thresh- 
old varied up and down by a factor of two from the default value. 

The impact of initial- and final-state radiation is estimated using a tt MC sample generated 
with POWHEG, instead of MadEvent/MadGraph. 

We estimate the effects of these systematic uncertainties on the expected t'F cross section limits 
by adding them to the limit calculation one at a time. The largest effects on the expected lim- 
its come from the normalizations of the EW background, the jet energy scale calibration, and 
the normalization of the tt background. All other uncertainties change the expected limits by 
insignificant amounts. In order to simplify the computational complexity of the limit computa- 
tion, we therefore consider only a limited set of systematic uncertainties in the limit calculation 
by assigning nuisance parameters to them: the integrated luminosity, normalization of the EW 
and tt backgrounds, lepton efficiency, jet energy scale, and parton-shower matching threshold. 
The additional effect of the other uncertainties is negligible. All of these except the lepton effi- 
ciency are treated as correlated in the combined result from the e+jets and the ^+jets channels. 

8 Results 

Figure |6] shows the observed and expected 95% CL upper limits on the t'F cross section for the 
e+jets (top), the f/+jets (middle) channels, and the combination of both channels (bottom). The 
95% CL lower limit for the t'-quark mass is given by the value at which the observed upper 
limit curve for the t'F cross section intersects the theoretical curve, also shown in Fig. |6] In the 
e+jets channel this happens for the 95% CL observed (expected) lower limits for a t'-quark mass 
of 490 (540) GeV. In the f/+jets channel the corresponding t'-quark mass limit is 560 (550) GeV. 
The combined observed (expected) limit from both channels is 570 (590) GeV. A comparable 
upper limit on the t' mass of 557 GeV was obtained recently by the CMS Collaboration using a 
dilepton channel [38|. 

9 Summary 

The results of a search for up-type fourth-generation quarks that are pair produced in pp inter- 
actions at v^s = 7TeV and decay exclusively to Wb have been presented. Events were selected 
in which one of the W bosons decays to leptons and the other to a quark-antiquark pair. The 
selection required an electron or a muon, significant missing transverse momentum, and at 
least four jets, of which at least one was identified as a b jet. A kinematic fit assuming t'F pro- 
duction was performed and for every event a candidate t'-quark mass and the sum over the 
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Figure 6: The observed (solid line with points) and expected (dotted line) 95% CL upper limits 
on the t't' production cross section as a function of the t'-quark mass for e+jets (top), ;^+jets 
(middle), and combined (bottom) channels. The ±1 and ±2 standard deviation ranges for 
the expected limits are shown by the bands. The theoretical t'f' cross section is shown by the 
continuous line without points. 
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transverse momenta of all decay products of the t't' system were reconstructed. No significant 
deviations from the Standard Model expectations have been found in these two-dimensional 
distributions, and upper limits have been set on the production cross section of such t' quarks 
as a function of their mass. By comparing with the predicted cross section for t't' production, 
the strong pair production of t' quarks is excluded at 95% CL for masses below 570 GeV under 
the model assumptions used in this analysis. This result and the one from f3E\ are the most 
restrictive yet found and raise the lower limit on the mass of a t' quark to a region where per- 
turbative calculations for the weak interactions start to fail and nonperturbative effects become 
significant. 
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